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Abstract 


This article presents a novel numerical approach to the solution of the 
nonlinear Kawahara equation. The desired approximations are obtained 
from the combination of Dickson polynomials and Taylor’s expansion. The 


combined approach is based on Taylor’s expansion for discretizing the time 
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derivative and Dickson polynomials for space derivatives. The problem 
will be converted into a system of linear algebraic equations for each time 
step via some suitable collocation points. Error estimation is presented 
after obtaining the approximate solution. The newly proposed technique 
is compared with some existing numerical methods to show the method’s 
applicability, accuracy, and efficacy. Two problems are solved to demon- 
strate the method’s power and effect, and the results are presented as a 


table and graphics. 


AMS subject classifications (2020): 35C11; 35G20; 65N22. 


Keywords: Taylor’s expansion; Dickson polynomials; Kawahara equation; 


Discretization. 


1 Introduction 


Partial differential equations (PDEs), including linear and nonlinear PDEs, 
are used to describe many real-world phenomena [9, 25, 16]. We must solve 
differential equations to see the nature of the background of these phenomena. 
In general, we are better at dealing with linear problems than nonlinear prob- 
lems. However, in practice, we face an increasing number of nonlinear prob- 
lems. Nonlinear partial differential equations (NPDEs) have a significant role 
in science and engineering [27, 26]. In the last few decades, many researchers 
have solved numerous numerical approximations for solving NPDEs. The 
radial basis function partition of the unity method was used by Nikan and 
Avazzadeh [22] to solve the sine-Gordon system. Bashan et al. [5] proposed 
a numerical solution of the complex modified Korteweg-de Vries equation by 
the differential quadrature method (DQM). Also, Baghan et al. [6] proposed 
a finite difference method combined with DQM for numerical computation 
of the modified equal-width wave equation. Many numerical methods, such 
as finite difference, decomposition, homotopy perturbation, finite element, 
radial basis function, finite volume, and B-spline collocation methods, have 
been developed in recent years to solve these NPDEs [23, 29, 33, 30, 6]. 
Nonlinear evolution equations are employed in many engineering and math- 


ematics disciplines, including solid state physics, fluid mechanics, chemical 
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physics, plasma physics, geochemistry, and chemical kinematics[14, 4]. One 
of the well-known nonlinear evolution equations is the fifth-order Kawahara 
equation, which appears in the theories of shallow water waves having surface 
tension, magneto-acoustic waves in a plasma, and capillary gravity equations 
[8, 31]. 


Consider the nonlinear Kawahara equation [16] is of the form 


Ut + Utls + UZs + U55 = O, (s, t) € (0, 1] x (0, T| ’ (1) 


based on the boundary conditions 


72 105 4 
i 2 
u (0, t) 169 + 169 sech” ket, (2) 
_ 72 105 4 
u(1,t) = 765 + F@g sech k(1+ ct), (3) 
as well as the initial condition 
72 105 4 
=--1 4+ 4 
u(s, 0) 169 + 169 sech” ks, (4) 
where k = Wis and c= 3. 


Researchers have developed a wide range of approximate algorithms and 
analytical techniques over the last few decades to solve the Kawahara equa- 
tion. In scholarly literature, there are numerous effective methods for ap- 
proximating the solutions of nonlinear PDE. The homotopy analysis method 
was used by Abbasbandy [1] to solve the Kawahara equation. Turgut and 
Karakoc [2] proposed a numerical approach based on the collocation method 
to solve the modified Kawahara problem. To address the Kawahara problem 
numerically, Rasoulizadeh and Rashidinia [28] developed a local radial ba- 
sis function-finite dimension (RBF-FD) meshless approach. Haq and Uddin 
[13] looked at the RBFs approximation method for the Kawahara problem. 
Soltanalizadeh [32] used the differential transformation approach to solve the 
Kawahara equation numerically. Lahiji and Aziz [21] solved the Kawahara 
equation numerically using spectral methods. Karakog, Zeybek, and Ak [15] 
to compute numerical solutions of the Kawahara equation using the septic 
B-spline collocation method. Zara et al. [34] applied the kernel smooth- 


ing method to approximate the modified Kawahara equation numerically. 
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Several scholars have also developed different approaches for solving this 


equation. 


In this article, we describe a newly developed, extremely simple numeri- 
cal technique based on Dickson approximation. There are some studies, such 
as [20, 19, 18], that employed Dickson polynomials as a basis for numer- 
ically solving integral equations, but the present method is different from 
them. This study seeks to develop a newly combined approximation method 
that combines Dickson polynomial approximation with Taylor-series approx- 
imation to solve the nonlinear Kawahara equation. We discretize the time 
variables after first estimating the time-dependent term with Taylor’s series. 
The underlying model problem is then assumed to have a solution for the 
novel Dickson series expansion of the unknown function at each time step. 
Then, using an appropriate set of collocation points, all associated unknown 
functions are represented in the Dickson matrix form. Error analysis of the 
combined technique is also presented here. There are few works on solv- 
ing polynomial approximation since it is difficult to solve the NPDEs. So 
we combine Taylor’s scheme for time derivative and Dickson polynomial for 
space derivative to create a new technique to solve the Kawahara equation. 
The obtained results have been compared with the same obtained by some 
existing methods. It manifests that the present method is more accurate than 
other methods and the computational time taken by the present technique is 


much less (less than a second). 


The rest of the paper is organized as follows: Section “Introduction” 
highlights the importance of the Kawahara equation and its application to 
a certain field of study. It also mentions the main objectives or goals of the 
research, as well as any background information required for understanding 
the next parts. The “Preliminaries” Section emphasizes Taylor’s time dis- 
cretization scheme and the Dickson polynomials. The next part discusses the 
“Error analysis ” to quantify how closely the numerical solution approximates 
the exact solution. The section “Numerical methods for the Kawahara equa- 
tion” details the suggested numerical technique for approximating solutions 
to the Kawahara equation. It also explains the step-by-step procedure of 
the numerical method, possibly making use of the earlier discussed Dickson 


polynomials and Taylor’s time discretization scheme. The proposed method 
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is used to solve nonlinear Kawahara equations problems about their accu- 
racy and applicability in Section “Illustrative Examples.” The contributions 
of the work are summarized in the “Conclusion” section, which highlights 
the achievements of the proposed method and identifies possible directions 


for further study. 


2 Preliminaries 


Taylor’s time discretization scheme: 


First, the NPDE is attempted to be discretized with respect to the time 
variable; that is, the interval [0,7] can be divided into (M + 1) grid points 
as 

O=:to <t) =At<--+-<tmn = MAt=T, (5) 


with the uniform step size At = t,, —t,_1. To obtain a discretization scheme 
that is time accurate, we get the following from the Taylor series representa- 


tion for u” = u(s, ty): 


n n 
—u Atuf 


At 2 


+ O(At’). (6) 


Dickson polynomial: 


The Dickson polynomial was first introduced by Dickson [11] in 1896; after- 
ward, it was rediscovered by Brewer[10] in 1961. For integer n > 0 and a 
is a commutative ring R with identity. The first kind Dickson polynomial 


Dry, (s, a) of degree n with the parameter a € R is defined by 


3 


(—a)?s"—-?? , n>1,-co<t<o. 


m—p 


The first few Dickson polynomials are 
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1. Do (s,a) = 2, 
2. Di (s,a) =, 
3. Dz (s, a) = 8? — 2a, 
4. D3(s,a) = 8° — 38a, 
5. Da(s, a) = s* — 487a, 
6. Ds (s,a) = 8° — 5s°a + 5sa?. 


They may also be generated by the recurrence relation for n > 2 as 
Dy (8,@) = 8Dn_1 (8,a) — ADn_2 (8,0), n > 2. The polynomials u = 


Dry, (8, @) read the second-order differential equation: 


(s? — da)u” + su’ — n?u = 0, n=0,1,2,3,.... 


Function approximation: 


Since f(s) € L?(R) is a continuous function defined over [0,1], then the 


Dickson polynomials can be represented as 


f(s) = dL eiDils); s € (0, ]], (7) 


where D,,(s,a) = D,(s) with any arbitrary constant a. 
In this case, it is necessary to figure out the unknown coefficients, c;, 
i=0,1,...,N. The Dickson vector Dy(s) and the unknown vector Cy will 


be introduced in the following ways: 


T 
Dy(s) = [Do(s) Dis) -.. Da(s)] 
and 
T 
Cyn = [eo Cl tea , 
We can express (7) in a compact representation using these vectors as 


follows: 
f(s) = Cy Dw(s). (8) 
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Additionally, by introducing the matrix WV, 
2 0 0 ..0 
0 + (5) (—a)° 0 re 
ical] Oe” 3(c)(-a)? «0 
7 0 3(7)(—a)} 0 sal) ’ 
0 
N(* N N_(¥t1)\/_,)\¥-1 
ee cytes 
and the monomial vector 
T 
Sn(s) = [1s s?... s*| ; 
we shall express Dy(s) as follows: 
Dy(s) = WS,,(s). (9) 


For finding the operational matrix of the derivative, a straightforward calcu- 
lation shows that 


d d 


(10) 
= WrS;,(s), 
where 
010...00 
002 ...00 
A= : 0 (11) 
00...0 NO 
00; 0090 (N+1)x(N+1) 


3 Error analysis: 


Kiirkcii, Aslan, and Sezer [19] employed the residual function in the Banach 
space method to do convergence research on the Dickson polynomial solution 


for integro-differential equations. We can find the Dickson polynomial-based 
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error analysis for approximation for the nonlinear Kawahara equation. Fi- 
nally, using the Dickson and Taylor’s approximation, we derive the error 


analysis of approximation. 


Theorem 1. [17] Suppose that H is a Hilbert space and that Y is a closed 
subspace of H such that dim Y < oo and 1, y2,---,Yn is any basis for Y. Let 
x be an arbitrary element in H and let yo be the most accurate approximation 


2 G(2,Y1 Yay Yn 
of x € Y. Then ||x — yo||; = Spat where 


(x, 2) (x,y1) (x, ye) «++ (2, yn) 
GE Yi Poses. Va) = (ya, 2) (Yr, 41) (Yrs ¥2) +++ (Ya, Yn) 


(YnsZ) (Yrs Y1) (Yns Y2) *** (Yrs Yn) 


Theorem 2. Let u € C(([0, 1]x[0, 1]) and let Y = Span {Do, Di, Do,..., Dn}, 
where D,’s are Dickson polynomials. If uy = c’¢(zx) is the most accurate 
approximation to u for particular time step At, then the error bound is as 


follows: 


k 


al 
G(u Do dD, Bes Dy) 2 
uU— UW < O(At)? 4 ( 2? 
| nila S O(At) G(Do,Di,..-;Dn) 
Proof. Let u(x,t) be the exact solution and let u(x,t,) be an approximate 
solution after time discretization. For a particular time step At, from (6), we 


have 
l|u(a,t) — ula, tn)llp < O(At)?. 


Let uy (2, tn) = ec o(x) be the most accurate approximate solution of u(z, tn) 
using Dickson polynomials for a particular time step t = t,. The error bound 


is as follows: 


Ilu(x, t) — un (2, tr)Ilp = lula, t) — ula, tr) + ue, tn) — Un (2, tr) Ilo 
< |lu(a, t) — u(@, tn lle + [lu tn) — Un (2, tn )Ilo 


G(u, Do, Di,..., Dn) 2: 
G(D6; Dioxsegd dy) 


< O(At)? + ( 
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4 Numerical method for Kawahara equation 


Consider the nonlinear Kawahara equation defined in (1)—(4) and the differ- 


entiate equation (1) with respect to t, we get 


Ur, = (Uss)e — (Uas)e — (wets )e- (12) 
n+1 nT 
By replacing first-order derivatives u? ~ Cae in all occurrences, we 


may write the equation as 


uf, = (=) (Cia (yu u( 2 = #8) 
At 5s At 3s At : At 


Atuiy = (uggs) — uss) — (use — UBy) — (UPFY — uy — u(ugh? — us). 


Next, using the time discretized form, (13) is inserted into the right-hand 
side of (6) to get 


At (14) 
(ut — aug — (ul? — ut) 
Now, discretizing the time-dependent times of (1), we get 
Uy = Ub, — UZ, — U"UL. (15) 
Again, replace (15)in the left side of (14); that is, 
n n Th, Mh arth — ag i n n n 
oy — hy aha (EE) Maen = Gant — uy 
(uj — wh.) — (ul — ul )uy — u(ut*? — at] 
(16) 
Equation (16) can be simplified as 
1 1 1 1 -1 
mae ae n+1 apt n+1 = n+1 acs n+1 
lag t+ gusle | 5 ]uz + [5] us. +[ 5 | Ube a7) 


1 1 u” 


n | 


55s 9 3s t At’ 
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The unknown functions u"+!,u?+1, u%t!, and uZ{! can be approximated 


by Dickson polynomials defined in (17) as 
i s) = CL DNs) = CLWS, (8). (18) 


Similarly, 


(s) 
ugs? = CHV Sy ( ) 
ue’ = CWS. (s) 

(s) 


ter =O re) 


Now, (17) can be described as follows: 


A3(s)CEWX*S,(s) + Ao(s)CE VAS, (s) 


(19) 
+ Ai(s)CZWAS,,(s) + Ao(s)CEVS,(s) = Fr(s), 
where 
al i 1, 1 
Ao(s) => At + gus? A;(s) =x7u , Ao(s) = Py 
1 en 1, u” 
A3(s) ~)5? Fy (s) —_ 9 bs 9 3s At’ 


and (19) can be simplified as 
CEW[A3(s)A° + Ao(s)3 + Ai(s)d + Ao(s)|Sn(s) = F,,(s). (20) 


Putting the collocation points s; = (+) in (20), we get a system of the 


linear equation as follows: 
Cy W[A3(si)A° + Ao(si)A? + Ai(s:)A + Ao(si)]Sn(si) = Fr(si) (21) 


fori = 0,1,2,...,N. Equation (21) is clearly a set of (N +1) linear equations 
in terms of (N + 1) unknown coefficients co,ci,¢2,...,cn that need to be 
found. 


Again, (21) can be expressed as a matrix equation as 


DynCy = Fh, (22) 
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where 
Dn := (®A3(A7)° + ®Aa(A7)3 + BA1(AT) + BAQ) UT 


and 


Am(so) 0 0 
0 Am(s1) 0 
A — 
0 0 vee Am(sn) (N+41)x (N41) 
fn(so) 
fn(S1) 
nel 
In(sy) (N41)x1 
From the given boundary condition, we have uy+1(0) = hor = 3 dhe 
103 sech* ket and unyi(1) = htt? = —B + 18 sech*k(1 + ct), using the 


matrix notation, it may be represented as 

D,(0)\Cw = hg", D,(0) = 8,,(30) 8", 

DADC rH he, Doll Sa 

Next, we substitute the first two rows of (22) by the boundary condi- 
tions [D,,(0);hgt"] and [D,,(1);h?**] for the convenience. From the above 


system, we can obtain the approximate value of the unknown coefficients 


Co, C1, C2,--.,¢n and hence obtain the solution of (18). 


5 Illustrative examples 


To demonstrate the performance of the Combined Taylor—Dickson Approx- 
imation (CTDA) technique, numerical simulations based on the nonlinear 
Kawahara equation for specific initial and boundary values are provided. 
We utilize an ASUS laptop with Processor (CPU) Ryzen 5 (3500U) and 
RAM(8GB) processing capabilities for conducting the computation. The 
MATHEMATICA software is used for the problem implementation and 3D 


Iran. J. Numer. Anal. Optim., Vol. 14, No. 1, 2024, pp 20-43 


31 


A new numerical approach to the solution of the nonlinear Kawahara... 


plotting of exact solution and approximate solution graphs with various val- 
ues of s, t, and At. 


Example 1. Consider the nonlinear Kawahara equation, (1)—(4). The exact 
solution of this problem is given byKaya and Al-Khaled [16] as u(s,t) = 
— 72 + 19 sech* k(s + ct), where k = Wi and c = 3%. 
At = 0.1, At = 0.01, At = 0.001, and N = 6, the absolute errors obtained by 


the present method for 0 < s < L = 1 is given in Table 1. Despite the fact 


Using time steps 


that there are few approximation methods for solving the Kawahara problem, 
the authors have compared their findings to the exact answer that has been 


verified. 


The maximum absolute errors solved by the RBF method [13] with the 
nonlinear Kawahara equation are compared in Table 2. Table 3 displays the 
absolute error for Example 1 determined using the CTDA technique for vari- 
ous values of the Dickson parameter a = 1,2, and the solution is not affected 
for different values of a. Figures 1 and 2 show the behavior of the exact 
and approximation solutions for N = 6 and At = 0.1, respectively. Figure 3 
shows a two-dimensional graphical representation of the exact and approx- 
imate solutions for t = 1 and At = 0.1. The graph represents how these 
possibilities compare visually. Figure 4 displays an “Absolute error graph” 
that shows the differences between the exact and approximate solutions for 
N = 6 and At = 0.1. In the literature, the computation time of the RBF 
method has not been provided. However, the proposed method provides 
more accurate results than the RBF method, and the computational time of 
the present method is much less than 0.267 seconds for N = 6. 

From Theorem 2, the maximum absolute error for N = 6 and At = 0.1 


is computed as 


dik 


= (0.1)? + 2.01111 x 1077 


O(At)? ae (S32) 


G(Do, D1,.--,Dn) 
= 1.00002 x 1077. 


Similarly, for At = 0.01 and At = 0.001, the maximum error can be calcu- 
lated as 1.00201 x 10-4 and 1.20111 x 10~®, respectively. So, this present 


method satisfies the error-bound condition for different time steps. 
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Table 1: Absolute error for Example 1 for N = 6. 


t | s At =0.1 At = 0.01 At = 0.001 
0.2 | 0.2 | 2.02289 x 10-® | 2.35173 x 10-8 | 2.66443 x 10-" 
0.4 | 2.2706 x 10-® | 1.01181 x 10-7 | 2.66443 x 10-10 
0.6 | 1.33716 x 10~® | 2.24034 x 10-8 | 2.38497 x 10710 
0.8 | 9.9872 x 10-® | 4.80536 x 10-8 | 2.66443 x 10-10 
0.4 | 0.2 | 3.82575 x 10-® | 4.60788 x 10-7 | 3.20035 x 10710 
0.4 | 5.57816 x 10-§ | 4.54266 x 10-7 | 2.24848 x 10-9 
0.6 | 1.9494 x 10-® | 1.72914 x 10-8 | 1.91345 x 10-9 
0.8 | 8.80955 x 10-® | 4.29541 x 10-8 | 1.91345 x 107-9 
0.6 | 0.2 | 4.51829 x 10-® | 5.54659 x 10-7 | 1.83046 x 10-9 
0.4 | 7.4988 x 10-® | 7.01482 x 10-7 | 1.6928 x 10-9 
0.6 | 5.48562 x 10- | 7.81906 x 10-7 | 1.29511 x 10~° 
0.8 | 3.23781 x 10-® | 2.91407 x 10-8 | 1.91345 x 10-9 
0.8 | 0.2 | 3.41209 x 10-® | 4.23962 x 10-8 | 8.60722 x 10-9 
0.4 | 6.12709 x 10-® | 6.07458 x 10-7 | 1.10209 x 10-10 
0.6 | 5.48562 x 10-® | 1.03009 x 10-7 | 8.43569 x 10710 
0.8 | 1.20464 x 10-® | 1.47048 x 10-8 | 1.91345 x 10-9 


32 


Table 2: Comparison of Loo errors of Example 1 for At = 0.1 and N = 6,8, 10. 


9.9872 x 10-4 
2.64883 x 10~® 
2.223943 x 107-7 

9.38 x 10-6 


Present method(N = 6) 

Present method(N = 8) 

Present method(N = 10) 
RBF(13] method(N = 51) 
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Table 3: Absolute error of Example 1 for Dickson parameter a = 1, 2. 
t x a=1 a=2 
At = 0.1 At = 0.01 At = 0.1 At = 0.01 
0.1 | 0.1 | 3.967 x 10-® | 7.253 x 10-° | 3.967 x 10-® | 7.253 x 10-8 
0.2 | 9.746 x 10-® | 1.098 x 10-8 | 9.746 x 10-® | 1.098 x 1078 
0.3 | 1.114 x 10-® | 9.507 x 1078 | 1.114 x 107 | 9.507 x 1078 
0.2 | 0.1 | 7.927 x 10-© | 1.463 x 10-8 | 7.927 x 10~-© | 1.463 x 1078 
0.2 | 2.022 x 10-® | 2.351 x 1078 | 2.022 x 10-® | 2.351 x 1078 
0.3 | 2.559 x 10-® | 2.382 x 1078 | 2.559 x 10-® | 2.382 x 1078 
0.3) 0.1 | 1.114 x 107° | 2.134 x 1078 | 1.114 x 107® | 2.134 x 1078 
0.2 | 3.012 x 10-® | 3.575 x 1078 | 3.012 x 10-® | 3.575 x 1078 
0.3 | 4.006 x 10~® | 3.966 x 10-8 | 4.006 x 10~® | 3.966 x 10-8 


Figure 1: Exact solution of Example 1 for At = 0.1 and N = 6. 
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Figure 2: Approximate solution of Example 1 for At = 0.1 and N = 6. 


0.195 - 


0.190 F 


0.185 £ 4 


u(s,t) 


0.180 | 


0.1755 | 


0.0 0.2 0.4 0.6 0.8 1.0 


Figure 3: Two-dimensional graph of Example 1 for At = 0.1 and N = 6 and t = 1. 
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| 0.00015 


} Absolute Error 
| 0.00010 


0.00005 


»} 0.00000 


Figure 4: Absolute error of Example 1 for At = 0.1 and N = 6. 


Example 2. Here, we consider a Kawahara-type equation given as 

Ut + UUs + UZs + Us — Us5s = 0. 

Whose exact solution is given as 

u(s,t) = —702 sech* k(s — 2 — ct), where k = ONaE and c = 72°. Using 
time steps At = 0.1, At = 0.01, At = 0.001 and the Kawahara equation with 
N = 8, the approximation for 0 < s < EL = 1 is given in Table 4. Despite 
the fact that there are few approximation methods for solving the Kawahara 
problem, the authors compare their findings to the exact answer that has 


been verified. 


The maximum absolute errors solved by the RBF method [13] with the 
nonlinear Kawahara equation are compared in Table 5. Table 6 displays the 
absolute error for Example 2 determined using the CTDA technique for vari- 
ous values of the Dickson parameter a = 1,2, and the solution is not affected 
for different values of a. Figures 5 and 6 show the behavior of the exact 
and approximation solutions for N = 8 and At = 0.1, respectively. Figure 7 
shows a two-dimensional graphical representation of the exact and approx- 
imate solutions for t = 1 and At = 0.1. The graph represents how these 
possibilities compare visually. Figure 8 displays an “Absolute error graph” 
that shows the differences between the exact and approximate solutions for 
N =6, At = 0.1. In the literature, the computation time of the RBF method 


has not been provided. However, the proposed method provides more accu- 
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Table 4: Absolute error for Example 2 for N = 8. 


t s At =0.1 At = 0.01 At = 0.001 
0.2 | 0.2 | 1.894 x 10-6 1.73465 x 1078 | 3.13242 x 1071! 
0.4 | 9.10089 x 10-® | 1.21102 x 10-” | 1.23812 x 1071° 
0.6 | 1.71853 x 10-® | 2.04784 x 10-8 | 2.06468 x 101° 
0.8 | 2.99047 x 10-® | 2.01218 x 10-8 | 2.01002 x 10~1° 
0.4 | 0.2 | 2.64883 x 10-© | 4.13838 x 107-7 | 4.36423 x 1071° 
0.4 | 1.32887 x 10-6 | 1.77252 %10-" | 1.81323 10-? 
0.6 | 2.59029 x 10-® | 3.0947 x 107-8 3.12158 x 10-9 
0.8 | 2.89771 x 10-© | 3.16247 x 1078 3.161 x 1079 
0.6 | 0.2 | 2.46179 x 10-® | 3.81904 x 10-7 | 4.03194 x 10-° 
0.4 | 1.28578 x 10-® | 1.71872 x 10-” | 1.75922 x 10-° 
0.6 | 2.5832 x 10-® | 3.09405 x 10-® | 3.12232 x 10-9 
0.8 | 2.99047 x 10-© | 3.27557 x 10-8 | 3.27585 x 107° 
0.8 | 0.2 | 1.5832 x 107° | 2.32987 x 107® | 2.46238 x 10-9 
0.4 | 8.23146 x 10-® | 1.110252 x 10-7 | 1.112919 x 10719 
0.6 | 1.70397 x 10-® | 2.04625 x 10-" | 2.0659 x 1071° 
0.8 | 2.03731 x 10-® | 2.223943 x 10-7 | 2.224081 x 10719 
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rate results than the RBF method and the computational time of the present 
method is much less that is 0.326 seconds for N = 8. 


6 Conclusion 


The Kawahara equation was solved using a combination of Taylor’s and Dick- 
son polynomial approximation. The authors employed derivative operational 
matrices as a tool to implement their chosen approximation methods. These 
matrices help in the numerical computation of function derivatives, which is 
necessary when dealing with differential equations. For small values of N, 
the current method yields very precise results with exact solutions. We get a 
more accurate solution if we increase the value of N. Also, taking small val- 
ues of At yields a more accurate result. The validity and convenience of the 


existing techniques were demonstrated using numerical examples. These ex- 
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Table 5: Comparison of Loo errors of Example 2 for At = 0.1 for N = 6,8, 10. 


Present method(N = 6) 


1.3288 x 1074 


Present method(N = 8) 


1.28578 x 10~® 


Present method(N = 10) 


4.36423 x 10-8 


RBF[13] method(N = 51) 


6.159 x 107° 


Table 6: Absolute error of Example 2 for Dickson parameter a = 1, 2. 
t x a=1 a=2 
At = 0.1 At = 0.01 At = 0.1 At = 0.01 
0.1] 0.1 | 1.118 x 107-© | 2.285 x 1078 | 1.118 x 107© | 2.285 x 1078 
0.2 | 3.220 x 10-® | 4.155 x 1078 | 3.220 x 10-® | 4.155 x 1078 
0.3 | 5.039 x 10-8 | 4.523 x 10-® | 5.039 x 10-® | 4.523 x 10-8 
0.2 | 0.1 | 2.188% 10-8 | 4.117 10-8 | 2.138 x 10-% | 4.117 x 10-8 
0.2 | 5.860 x 10~® | 7.590 x 10-8 | 5.860 x 10~® | 7.590 x 10-8 
0.3 | 8.355% 10-° | 9.451 « 10-* | 8.355 x 10-° | 9.451 « 10-8 
0.3. | 0.1 | 2.835 x 10-° | 5.463 x 10- | 2.835 x 10-® | 5.463 x 10-8 
0.2 | 7.847 x 10-® | 1.019 x 10-8 | 7.847 x 10-® | 1.019 x 1078 
0.3 | 1.133 x 107® | 1.296 x 1078 | 1.133 x 107® | 1.296 x 1078 


Figure 5: Exact solution of Example 2 for At = 0.1 and N = 8. 
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Figure 6: Approximate solution of Example 2 for At = 0.1 and N = 8. 


exact 


0.58 b —--— aprox 


u(s,t) 


fen a a rr 


0.0 0.2 0.4 0.6 0.8 i 


Figure 7: Two-dimensional graph of Example 2 for At = 0.1, N = 8, and t=1. 
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Figure 8: Absolute error of Example 2 for At = 0.1 and N = 6. 


amples also proved the accuracy and efficiency of the proposed method. This 
method established a balance between accuracy and computational efficiency, 
which is important when dealing with challenging mathematical problems. 
This method can also be used to solve various NPDEs derived from physical 
models, which are subject to specified initial and boundary conditions. The 
method’s future direction is that it can be used for NPDE with higher-order 


time derivatives. 
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